Localization of spin mixing dynamics in a spin-1 Bose-Einstein condensate 



O 
(N 



m 
(N 



Wenxian Zhang, ^ Bo Sun,^ M. S. Chapman,'^ and L. You'^ 

^ The Key Laboratory for Advanced Materials and Devices, 
Department of Optical Science and Engineering, Fudan University, Shanghai 200433, China 
^Department of Physics, Auburn University, Auburn, Alabama 36849, USA 
'^School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430, USA 

(Dated: March 24, 2010) 

We propose to locahze spin mixing dynamics in a spin-1 Bose-Einstein condensate by a tem- 
poral modulation of spin exchange interaction, which is tunable with optical Feshbach resonance. 
Adopting techniques from coherent control, we demonstrate the localization/freezing of spin mixing 
dynamics, and the suppression of the intrinsic dynamic instability and spontaneous spin domain for- 
mation in a ferromagnetically interacting condensate of *^Rb atoms. This work points to a promising 
scheme for investigating the weak magnetic spin dipole interaction, which is usually masked by the 
more dominant spin exchange interaction. 
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Dynamic localization is ubiquitous in nonlinear sys- 
tems, both for classical dynamics as in an inverted pen- 
dulum with a rapidly modulating pivot [ij or an ion in 
a Paul trap 0, and for quantum dynamics like a one- 
or two-dimensional soliton in a Bose-Einstein condensate 
(BEC) when the attractive interaction strength is rapidly 
modulated [^-Q. It is often used to stabilize a dynami- 
cally unstable system. 

Spin mixing dynamics of a spin-1 atomic condensate 
are dynamically unstable [3| when the spin exchange in- 
teraction is ferromagnetic, i.e., favoring a ground state 
with all atomic spins aligned. When confined spatially, 
the unstable dynamics is known to cause formation of 
spin domain structures @, [SS] . For many applications of 
spinor condensates jfrom quantum simulation to preci- 
sion measurement [lO|, it is desirable that spin domain 
formation is suppressed. In addition, atomic spin dipo- 
lar interactions, although weak when compared to typical 
spin exchange interactions, induce intricate spin textures 
that are difficult to probe when masked by spin domain 
structures. Thus the suppressing/freezing of the unde- 
sirable dynamics from spin exchange interaction is also 
important for investigating the effect of dipolar interac- 
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Compared to conventional magnets in solid states, a 
spin-1 BEC has one unsurpassed advantage: its spin ex- 
change interaction between individual atoms can be pre- 
cisely tuned through optical (as well as magnetic) Fes- 
hbach resonances 13l4l9l|. By adjusting the two s-wave 
scattering lengths ag and 02 of two colliding spin-1 atoms 
via optical means, the spin exchange interaction strength, 
characterized by C2 = 47rft^(a2— ao)/3Af with M the mass 
of the atom, is tunable. Analogous to an inverted rigid 
pendulum with a rapidly oscillating pivot, a fast tem- 
poral modulation of the spin exchange interaction can 
localize the spin mixing dynamics, equivalent to a sup- 
pressing/nulling of the spin exchange interaction. 

This study is devoted to a theoretical investigation of 



spin dynamics in a spin-1 BEC under the temporal mod- 
ulation of the spin exchange interaction. As an applica- 
tion, we illustrate the suppression of the dynamic insta- 
bility and the resulting prevention of spin domain forma- 
tion in a condensate with ferromagnetic interaction. The 
proposed scheme to control the spin exchange interaction 
will potentially provide a substantial improvement to the 
accuracy of several envisaged magnetometer setups and 
to enable cleaner detections of dipolar effects. 

For both spin-1 atoms ^^Rb and ^'^Na, popular ex- 
perimental choices, their spin-independent interaction 
strength, characterized by cq = 47r?i^(2a2 + ao)/3M, is 
two to three orders of magnitude larger than \c2\ |20l - 
22| . This ensures the validity of single spatial mode ap- 
proximation (SMA) (isl - ETt when the number of atoms 
is small and the magnetic field is low. The spin degrees 
of freedom and the spatial degrees of freedom become 
separated within the SMA. This allows one focus on the 
most interesting spin dynamics free from density depen- 
dent interactions. 

Within the mean field framework, the spin dynamics of 
a spin-1 condensate under the SMA is described by 
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where pi (i = -|-, 0, — ) is the fractional population of com- 
ponent \i), (X^i Pi — ITT- = P+ — P- is the magnetiza- 
tion in a spin-1 Bose condensate, a conserved quantity. 
9 is the relative phase fisf . (t>{r) is a unit normalized 
spatial mode function under the SMA determined from a 
scalar Gross-Pitaevskii equation with an s-wave scatter- 
ing length of a2. As before, the effective spin exchange 
interaction is given by c{t) = C2{t)N J df\<j){f)\'^ , albeit 
the time dependence, with N the total number of trapped 
atoms. Although the system dynamics ([T]) does not con- 



serve the total spin energy 




0.7 


\fW\il 




£{t) ^ c{t)pQ 


(1 - po) + - Po)^ - Tm? cos e 


,(2) 


c? 0.5 


1 n i ! ' 
* ,/ 'J li 
\;l M'l 


' ! 1 \ 1 

/ Wi' i 


due to the temporal modulation, the transversal spin 
squared fj_ — fx+ fy = 2£{t)/c{t) remains conserved 


0.3 


Is,' ■ W 1 : 

\ 1 \ I ■ 


\\ i ^ / 



and is determined solely by the initial condition. Because 
£{t) and c(t) are modulated exactly in the same manner, 
replacing 9 with equation ([1]) is further simplified to 
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fiPu ~ Po)(PO - Pd), 



(3) 



where p^A = fl [l ± V^~P) IV with p = fl+m\ 
Pu(d) takes the +(— ) sign, denoting the largest (smallest) 
value of Po along the orbit and satisfies polpu d ~ 0- It is 
straightforward to find the solution 



Po{t) 



Pu+ Pd , Pu - Pd 



■ sm[7 - 



P{s)ds], (4) 



where P{t) = ±2c{t)f /h is the frequency of the periodic 
spin evolution and the ± sign denotes the forward and 
backward evolutions, respectively, and 
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Po{0)-[{pu + Pd)/2] 

ViPn - Po(0)][po(0) - Pd] 



is given by the initial values of po and 0. The above 
solution (j4]) is valid for an arbitrary temporal modulation 
function c{t). 

To control the spin dynamics, we consider several sim- 
ple but practical modulations in the following. Based 
on these examples, we demonstrate that spin dynamics 
with a modulated c is very different from the free dy- 
namics and understand how a temporal modulated c{t) 
affects the spin dynamics. 

First we consider a sinusoidal modulation with c{t) = 
dcos{nt). d and are respectively the modulation am- 
plitude and frequency. The solution Eq. Q then becomes 



Po(t) 



Pu + Pd ^ Pu - Pd . r , . 

1 sm[7 + 7/sm(i2<)J, (5) 



with rj = zL2df/hn. The corresponding results are il- 
lustrated in Fig. [U for n = 0, 1/2, 1, and 2. The case 
of f2 = is simply the free evolution without modula- 
tion with a period To = 27r//3 determined by the initial 
condition For other cases, irrespective of the values 
for r2, the frequency of oscillation is always ft and the 
corresponding period is 2Tr/n,. As shown in Fig. [TJ the 
oscillation amplitudes show two distinctive regions: one 
for r2 < = 2df /nh (i.e. rj > n) where the amplitude is 
the same with/ without modulation; another for D, > flc 
(i.e. T] < tt) where the amplitude decreases with mod- 
ulation frequency 51. In this latter region, it is easy to 
check that A — (p„ — Pd)(l — cos ?/)/2 for the case shown 
in Fig. m 
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FIG. 1: (Color online) Left panel: Time dependent fractional 
population poit). The black dotted line refers to the free 
evolution without modulation, while the blue dash-dotted, 
red dashed, and blue solid lines are for O = 1/2, 1, and 2, 
respectively. The parameters and initial conditions are h= 1, 
d = -1, m = 0, po(0) = 0.7, 61(0) = 0, and 0.^ ~ 0.58. Right 
panel: The dependence of oscillation amplitude A of po on 
the modulation frequency fi. 




FIG. 2: (a) Orbits without modulation; Orbits with modu- 
lation for (b) Q, = 1/2, (c) = 1, and (d) Q. = 2. Parameters 
are the same as in Fig. [1] 



Figure [2] illustrates the orbits for the corresponding 
spin dynamics. A full orbit is occupied if 17 < i7c but 
only partial orbits are occupied when $7 > Jlc- The occu- 
pied portion deceases when f2 increases. Spin dynamics 
for the first half period is reversed during the second half 
period evolution, irrespective of the values for f2 ^ 0. 
This reversal is responsible for a more robust modulated 
dynamics against various noises as noises are not reversed 
and their effect can be averaged zero according to coher- 
ent control theory [30;]. 

Next we consider a periodic square function modula- 
tion with 



c{t)^ 



d, n{w + t) <t < n{w + r) + w, 
0, n{w + t) +w <t < {n + l){w + r), 



(6) 



for n = 0,1,2,---. The spin dynamics is halted com- 
pletely if n{w + r) + w < t < (n + l){w + t) and is 
unmodulated if n(w + r) < t < n{w + t) + w. The corre- 
sponding plots for Po and 6 display interesting step like 
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FIG. 3: (Color online) (a) Time-dependent fractional pop- 
ulation po and (b) time-dependent relative phase 9 for the 
periodic square modulation Eq. ([6]) with t — w — 1/4 and 
d= -1. 



features as shown in Figs.[3lja) and[3l^b). 

The modulation dynamics considered above offers 
many interesting possibilities. A direct application is to 
remove a dynamical instability observed in a ferromag- 
netically interacting spin-1 condensate 0-S 31, 32|. This 
instability is removed whenever the imaginary part of the 
eigenfrequency for the corresponding Bogoliubov excita- 
tion becomes zero in a modulation cycle. We show this 
instability is indeed suppressed in the following for the 
cosine modulation C2(t) — dcos(flt). This suppression 
inhibits the spontaneous formation of spin domains. 

We now consider our system of a homogeneous *^Rb 
spin-1 condensate starting from an off-equilibrium initial 
state Q- The averaged spin / = {Fx)x + {Fy)y -\- {Fz)z 
is conserved where F^^y^z are spin-1 matrices. Starting 



from any stationary point, the evolution of the collective 
excitations takes a compact form 



(7) 



where x = {5^ +,6'^{),5'^ ^,6'^\,5'lil,5'li*_Y with devi- 
ations 5"^ j and 5'^* from the stationary solution 7]. 

The general solution to Eq. © is x{t) = U{t, 0)f (0) 
where U{t2,ti) = T exp[-(i/^) J^^ dsMis)] with T the 
time ordering operator. For periodic modulation, the so- 
lution during t G [pT, (p -|- l)r] is further simplified to 
x{t) ^ U{t,pT){UTfx{0), where p = 0,1,2,- •• indexes 
the number of periods. T — 27r/f2 is the period of modu- 
lation, and Ut = U{T,0) = Texp[-{i/h) dsM{s)] is 
the evolution operator for a complete period. \x{t)\ will 
grow (decay) exponentially with p if the modulus of Ut 
are larger (smaller) than unity. Diagonalizing Ut and 
rewriting it as Ut = V'' exp{—iTF)V, where the Fluo- 
quet operator F is a 6-by-6 diagonal matrix, the crite- 
ria for stable dynamics reduces to a vanishing imaginary 
part of Fq {q = 1,2, ••• ,6). Unstable dynamics arises 
if the imaginary part is not zero, while stable dynamics 
emerges if the imaginary part of all Fq is exactly zero. 

The inset of Fig. |4] shows the imaginary part of a typ- 
ical spectrum for the system under a cosine modulation 
of C2. We focus on the most unstable mode which in 
principle dominates the unstable dynamics. Compared 
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FIG. 4: (Color online) The dependence of the wave vector 
k- for the most unstable mode on the modulation frequency 
fl. The inset illustrates the imaginary part of a typical Bo- 
goliubov spectrum under modulation (blue double contours). 
For comparison, black dashed lines denote the results without 
modulation. 



to the modulation free results (in dashed lines), we find 
the most unstable mode (in double contours) is not only 
suppressed in amplitude but also shifted to larger wave 
vector. The dependence of fc_ on the modulation fre- 
quency is illustrated in Fig. |4l The almost indepen- 
dence of A;_ on Q at small values contrasts with a strong 
monotonic increase at large values of il. 

The emergence of spin domains is prohibited due to the 
modulation. On one hand, the suppression of F implies a 
smaller effective spin exchange interaction thus a longer 
effective spin healing length ^; On the other hand, the up- 
shift of fc_ means shorter wave length A = 27r/fc_. If the 
unstable mode wave length (potentially domain width) 
is smaller than the spin healing length, the condensate is 
able to heal by itself. The domain structure would never 
appear if ^ exceeds A. 

Because of the modulation, however, the resulting dy- 
namics becomes completely different from the case of 
C2 = or no spin dynamics at all. The modulation does 
not stop spin mixing dynamics, i.e., as we continue to 
observe spin waves which is nominally disguised in exper- 
iments by the spontaneously formed spin domains fs, ■ 
Furthermore, we expect the modulated spin dynamics 
to be robust against various experimental noises because 
the periodic modulation effectively cancels uncorrelated 
noises from alternating modulation periods. 

Finally we confirm the above conclusions for a trapped 
spin-1 condensate with full numerical sirnulations. We 
adopt experimental parameters as in Ref. [7| : The initial 
conditions are as in the experiment Q, with ®^Rb con- 
densates [po{0) = 0.744, 61(0) = 0, TV = 2.0 X 10^ and 
TO = 0], in a trap Ve^t{r) = {M/2){ujIx^ + w^y^ + a;2^2) 
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FIG. 5: (Color online) (a) Localization of the spin dynam- 
ics for a trapped spin-1 condensate in 3 cases: (i) mod- 
ulation starts at t = (red dashed line); (ii) no modu- 
lation (black dash-dotted line); and (iii) modulation starts 
at t = 6 (blue solid line). Spatial distribution m{z) — 



above three cases: (b) 
case (ii). a 



|$_(r, «)! ) at different times for the 
case (i); (c) — case (iii); (d) — 
Trivial and fiat m{z) at early times 



{t < 250) has been omitted. Dynamical instability induced 
spontaneous domain formation is prohibited by the modu- 
lated spin exchange interaction. 



with uj^ ^ ujy = (27r)240 Hz and = (27r)24 Hz. The 
modulation function is C2(t) — c? cos (fit) with d = C2, 
and — ujz which is about 3 times larger than the free 
spin evolution frequency 2tt/Tq. Two cases will be con- 
sidered: (1) The modulation is applied immediately (po 
oscillates around /5„); (2) The modulation is turned on 
at i = 6 (l/ujz) (po oscillates around pd). 

The results from numerical simulations are shown in 
Fig. [S] With modulation, the spin dynamics is clearly 
localized as (same for p^) oscillates with a smaller 
amplitude around its initial value, in contrast to the large 
amplitude unmodulated result [panel (a)]. In addition, 
the unmodulated dynamics shows domain structures af- 
ter t « 300, due to the intrinsic dynamical instability. 
While for both modulated cases, no spin domains are 
observed. Thus temporal modulation of spin exchange 
interaction does suppress the intrinsic dynamical insta- 
bility and prohibit spontaneous domain formation. In our 
extensive numerical simulations, we find that when addi- 
tional white noises are added, spin domains arc found to 
arise quicker for the unmodulated case; yet for the two 
cases with modulations, almost the same behaviors are 
observed as if no noise were added. 

Although the life time of the condensate at opti- 
cal Feshbach resonance is reduced dramatically in *^Rb 
gases 15|, we notice that there exists at least one magic 
window of relative frequency, e.g., between resonance /3 



and 7 in Fig. 7 of Ref. [l5|, where the condensate life 
time lasts several hundred milliseconds and C2 changes 
sign. On the other hand, the spin domain emerges in a 
time scale typically shorter than 100 ms 33j. Therefore 
the suppression of the domain formation in ®^Rb conden- 
sate is experimentally feasible. 

In summary, we propose to localize the spin mixing 
dynamics in a spin-1 condensate by temporally modulat- 
ing the spin exchange interaction. For condensed atoms, 
the modulation can be facilitated with the technique of 
optical Feshbach resonance [H, ■ We demonstrate the 
suppression of the intrinsic instability thus the inhibition 
of spontaneous spin domain formation in a fcrromagnet- 
ically interacting spin-1 Bosc condensate, such as *^Rb 
condensate in the _F = 1 manifold. In addition, the ef- 
fective freezing of spin mixing dynamics due to spin ex- 
change interaction provides a cleaner approach to inves- 
tigate magnetic spin dipolar interaction effect in a *^Rb 
Bose condensate [ll|, [l^ . 
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